Effect of quantum confinement of surface electrons on an atomic motion on 

nanoislands 
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Ab initio calculations and kinetic Monte Carlo simulations demonstrate that the quantum confine- 
ment of surface electrons to nanoislands can significantly affect the growth process at low tempera- 
tures. Formation of empty zones and orbits of an adatom motion is demonstrated for Cu nanoislands 
on Cu(lll). 



Noble (111) metal surfaces support Shockley surface 
states^. These electronic states arise in the inverted 
L gap of a metal band structure. Electrons occupying 
the Shockley-type surface states form a 2D nearly free 
electron gas confined in the vicinity of the top layer—. 
Surface electrons scatter at point defects^, adsorbates^, 
steps^, leading to standing waves of local density of states 
(LDOS), which can be detected using scanning tunnel- 
ing microscope (STM) technique. These standing-wave 
patterns contain important information about properties 
of scattering sites and interaction between the surface 
electrons and scattering centers. Particularly, scattering 
of the surface states at adatoms leads to the standing- 
wave patterns around the adsorbates and to an indirect 
long-range interaction (LRI) between adatoms^. While 
the first experimental evidence of such type of interac- 
tions has been reported 35 years ago-, quantitatively 
they have been resolved only recently by means of low- 
temperature STM studies^i^. These remarkable obser- 
vations have opened a door to creation of new types of 
macroscopic-ordered 2DiS and IDii nanostructures sta- 
bilized by the surface electrons. 

Several fascinating phenomena occur if surface elec- 
trons are confined to closed nanostructures, like corrals, 
vacancy holes or clusters. Experiments of Manoharan 
et al.fi^ and later ah initio calculations^ have demon- 
strated that the quantum confinement inside corrals in- 
duces a mirage effect. It is also possible to tailor spin- 
polarization of surface-state electrons and exchange in- 
teraction between magnetic adatoms within the con- 
fined nanostructures^^. Confined surface electrons inside 
quantum resonators alter diffusion at low temperatures, 
leading to atomic self-organization^^. To the best of our 
knowledge Li et al.fi^ performed the first quantitative 
investigation of the quantum confinement of surface elec- 
trons on nanoscale Ag islands on Ag(lll) by means of 
STM. These studies revealed the validity of the confine- 
ment picture down to the smallest of island sizes. It was 
proposed^^ that the quantum confinement can arise not 
only on Ag(lll), but on Cu(lll) and Au(lll). Observa- 
tion of quantized electronic states in vacancy islands on 
Cu(lll) has been recently reportedii. 

In this paper we demonstrate the effect of confined 
surface electrons on atomic motion in close proximity 



and on top of nanoislands. The quantum confinement is 
studied by means of the first principles Korringa-Kohn- 
Rostocker (KKR) Green's function method. Our kinetic 
Monte Carlo (kMC) simulations reveal that confinement- 
induced electronic states around and on top of nanois- 
lands significantly affect atomic diffusion. Formation of 
empty zones and orbits of adatom motion are shown. We 
demonstrate that the quantum confinement dramatically 
effects the growth process of nanoislands at low temper- 
atures. 

Within our study calculations of the electronic inter- 
action between an adatom and a nanoisland are per- 
formed using density functional theory in local spin den- 
sity approximation by means of the KKR Green's func- 
tion methodic. The basic idea of this method is a hier- 
archical scheme for the construction of the Green's func- 
tion of nanostructures on a surface by means of successive 
applications of Dyson's equation. The bulk, surface and 
impurity problems are consequently treated with a per- 
turbative approach. At each stage a fully self-consistent 
Green's function is obtained, which is then used as a ref- 
erence for the next step. We treat a surface as a 2D 
perturbation of an ideal bulk with a slab of vacuum. 
Taking into account the translational symmetry of the 
surface geometry, the Green's functions are formulated 
in momentum space. Adatom and nanoisland are con- 
sidered as a perturbation of a clean surface. These cal- 
culations are performed in real space. Previous studies 
have demonstrated that this method allows one to obtain 
with a good accuracy a spatial modulation of LDOS and 
interaction energy in closed nanostructure a ^'^i I'^d . 

As a model system we consider a hexagonal Cu nanois- 
land placed in fee hollow sites on a Cu(lll) surface 
and study diffusion of a Cu adatom near and on top of 
it. Each edge of the island has length of 11 atoms (2.6 
nm). Two different kinds of close-packed steps can be 
distinguished in the nanoisland: the {100}-microfaceted 
{A) step and the {lll}-microfaceted {B) step. We ne- 
glect the elastic interaction between the adatom and the 
island^^. 

At first, we concentrate on the quantum confinement 
on top of the island. Fig.l presents the 2D potential 
energy map of a Cu adatom placed in different hollow 
sites. The interaction potential exhibits an oscillatory 
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FIG. 1: (Color online) The 2D map of the interaction poten- 
tial between a Cu nanoisland and a Cu adatom placed on top 
of it. Semi-transparent white circles demonstrate the edges 
of the island and 1/3 part of it. The rows of hollow sites are 
numerated from the steps (numbers 1-4). 
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FIG. 2: (Color online) Atomic diffusion on top of the nanois- 
land: results of the kMC simulations for a single Cu adatom. 
The probability to find a randomly deposited Cu atom is cal- 
culated at 60 K. Semi-transparent white circles demonstrate 
the edges of the island and 1/3 part of its surface. The loga- 
rithmic scale is used. 



behavior if one recedes from an edge of the nanoisland 
towards its center along the perpendicular bisector to the 
edge. When the Cu adatom is located in the first row of 
hollow sites near an A or a, B step (Fig.l), the interac- 
tion is repulsive being larger than 30 meV. We attribute 
this effect due to redistribution of the electron-charge 
density at steps as was suggested a long time ago by 
Smoluchowski^^ . When the Cu adatom is located in the 
second row of hollow sites near B step the interaction is 
also repulsive being 10 meV. Our calculations reveal the 
attractive minima of -19 meV at 1.5 A from an A step 
(the second row of hollow sites) and -20 meV at 2.9 A 
from a B step (the third row). Difference in the positions 
of the minima with respect to the island edges appears 
due to the different types of atomic packing at A and B 
steps. The existence of attractive potential in a vicinity 
of a descending Cu step is confirmed by experimental ob- 
servations of Repp et al.jS. One of the most fascinating 
effects of confined surface electrons is quantum interfer- 
ence at the corners (black spots in Fig.l). The attractive 
potential in these spots is deeper than at steps being - 
25 meV. Similar phenomenon has been already observed 
for the scattering of surface-state electrons in a 'hand- 
made' triangular corral constructed of Ag adatoms on a 
Ag(lll) surface^^. 

In order to give a clear evidence that atomic motion 
on top of the nanoscale island is substantially modified 
compared to that on an open surface, we perform the 
kMC study2i of behavior of a Cu adatom randomly de- 
posited on top of the island. The hop rate of an adatom 
from site k to site j is calculated using the formulae: v — 
i/oea;p(— -Efc—tj/fc^T), where T is the substrate tempera- 



ture, is the attempt frequency (10^^ sec^^ (see Refi^)) 
and fcs is the Boltzmann factor. The infiuence of the 
electronic interaction 'adatom-island' (Fig.l) is included 
in the diffusion barrier: E^^^ = Ed + 0.5(£'j — Ek)^ 
where Ejj is the energy barrier for a Cu atom on a clean 
surface (40 meV^^) and £'fc(j) is the magnitude of inter- 
action, when the adatom is in the hollow site k{j). Pre- 
vious studies2^ have confirmed the applicability of this 
approach. 

The adatom probability distribution on top of the is- 
land is shown in Fig. 2. Contrary to the traditional view, 
the adatom does not diffuse as it does on a flat surface. 
Regions marked with the red color correspond to the 
energy minima of -25 meV in the corners of the island 
(Fig.l). The probability to find the adatom in these six 
spots is significantly higher than that one along the island 
edges. At the same time a localization of the adatom in 
a narrow stripe along the steps is more preferable than in 
the center of the island. Formation of empty zones and 
orbits of adatom motion takes place. 

Now we turn to the effect of the quantum confinement 
around the island. The 2D potential energy map caused 
by the substrate-mediated electronic interaction between 
the adatom and the nanoisland is presented in Fig. 3. The 
interaction has an oscillatory behavior, decaying to zero if 
one recedes from the island along the perpendicular bisec- 
tor to a step edge. The interaction is strongly attractive 
(about -1.4 eV), when the Cu adatom is at 2.2 A from 
the step; it corresponds to aggregation of the adatom and 
the island. Our calculations reveal a concentric repulsive 
ring of 25 meV at 4-7 A from the edge of the island. 
This barrier prevents atomic diffusion towards the steps 
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FIG. 3: (Color online) 2D interaction potential between a Cu 
nanoisland and a Cu adatom placed in its vicinity. The edge 
of the island is marked with the white circles. 



at low temperatures. There is a ring of attractive in- 
teraction of -5 meV at 8-14 A, and a ring of repulsive 
interaction of 2 meV at 15-20 A. The next minima and 
maxima of energy are less than 1 meV. Concentric rings 
of negative and positive energies (Fig. 3) lead to the for- 
mation of empty zones and 'allowed' orbits. Particularly, 
at low temperatures a migrating adatom prefers to be at 
about 9 A from the steps, while the probability to find 
an adatom at 4-7 A from the edge of the nanoisland is 
very low. 

The quantum confinement has a dramatic effect on 
growth of small clusters at low temperatures. To illus- 
trate this statement we perform the kMC simulations of 
self-organization of Cu atoms on Cu(lll) at 14 K for dif- 
ferent coverages of deposited atoms po- Within our study 
the LRI between Cu adatomsSS is described by a pair- 
wise summation. Previous studies have confirmed that 
this approximation is well-justified at large interatomic 
separations^^. 

In Fig. 4 we present histograms of statistically averaged 
density of adatoms on top of the island and around it 
as function of the distance between the adatom and the 
nearest step. Positive distances r correspond to the area 
outside the island, while negative ones to that on top of 
it. At coverage po = 1% of ML (Fig. 4a) there is only one 
peak of the density on top of the island. It indicates that 
all adatoms are localized at separations of 1.5-3 A from 
the steps, forming the first inner orbit. The first peak 
outside the island (at r = 2 A) corresponds to atoms 
nucleated to the island during deposition. The second 
maximum of atomic density outside the island is located 
at r = 9 A: these atoms form the first outer orbit. The 
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FIG. 4: Surface density of adatoms p(r) on top of the nanois- 
land and around it as a function of the distance r to the 
nearest island edge. Function p(r) is a number of adatoms 
over the square of a hexagonal concentric ring (marked with 
a shadow). Positive values of r correspond to the region out- 
side the island, while negative ones to that on top of it. The 
results of the kMC simulations of self-assembly of Cu adatoms 
at a presence of Cu nanoisland on Cu(lll) are demonstrated. 
The Cu adatoms are deposited simultaneously at 14 K; the 
coverage po is: a) 1% of ML, b) 3.5% of ML, c) 5% of ML. 
In Figs. b),c) peaks at r = -2 A are rescaled by factor 0.5. 
Dashed lines correspond to the uniform atomic density at a 
given coverage. 



probability to find an adatom on the second orbit (at r 
— 17-22 A) is much smaller. The electronic interaction 
'adatom-nanoisland' at such distances is less than 2 meV 
(Fig. 3), and the second orbit is stabilized by the LRI with 
the adatoms placed on the first orbit. 

Our simulations indicate that effect of the quantum 
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confinement on atomic self-organization remains signif- 
icant in a wide range of coverages. The results for po 
= 3.5% of ML are presented in Fig. 4b and for po = 5% 
of ML are given in Fig. 4c. With increasing the cover- 
age the probability of the occupation of the first orbit on 
top of the island (at A) increases. The second and 
the third inner orbits (at r~-12 A and at A re- 

spectively) are formed. The similar phenomenon occurs 
around the island: adatoms form orbits at A and 

r^28 A, while the atomic density between the orbits is 
low. Two Cu atoms located on the same orbit are sep- 
arated by 11-12 A, corresponding to the position of the 
local minimum of the LRI between two Cu adatoms^^. 

Finally, we would like to note that empty zones were 
found in several experimental works of Ehrlich and co- 
workers. In Ref.i-22.ilS diffusion of an Ir adatom around an 
Ir cluster on Ir(lll) was studied. Empty zone was re- 
solved at separations of 3-7 A from the edge of the clus- 
ter, while the preferable distance between the migrating 
atom and the step edges was 9-12 A. With increasing the 
temperature the adatom overcame the empty zone and 
nucleated with the cluster. The anisotropy diffusion was 
also observed for an Ir adatom placed on top of an Ir clus- 
ter on Ir(lll). The Ir atom was localized near the edge 
of the cluster with a higher probability than in its central 
region22i2^. Diffusion of a Pt atom on top of a Pt cluster 



on Pt(lll) was a subject of study in Refi^. It was found 
that empty zone (of width of 7 A) divides the region with 
uniform atomic distribution in the center of the cluster 
from the orbit near steps. During migration along the 
step edges, the preferable location of a Pt adatom was 
in the corners of the cluster. Results of our study for 
Cu suggest that phenomena observed by Ehrlich and co- 
workers could be promoted by the quantum interference 
of electrons scattered on the cluster, on one hand, and 
on the adatom, on another hand. However, in order to 
give a precise answer to this question, further theoretical 
investigations are required. 

In conclusion, we have performed ab initio study of 
the quantum confinement of surface electrons on top of 
nanoscale Cu islands on Cu(lll) and around them. For- 
mation of concentric rings with positive and negative in- 
teraction energies is revealed outside and on top of a 
nanoisland. This phenomenon substantially modify the 
atomic diffusion at low temperatures, leading to empty 
zones and orbits of adatom motion. Our results indicate 
on a profound role of the quantum confinement of sur- 
face electrons during growth of metal nanoclusters at low 
temperatures. 

This work was supported by Deutsche Forschungsge- 
meinschaft (SPP1153, SPP1165). 
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